Quantum Discord in a spin-1/2 transverse XY Chain Following a Quench 
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We report a study on the zero-temperature quantum discord as a measure of two-spin correlation 
of a transverse XY spin chain following a quench across a quantum critical point and investigate 
the behavior of mutual information, classical correlations and hence of discord in the final state 
as a function of the rate of quenching. We show that though discord vanishes in the limit of very 
slow as well as very fast quenching, it exhibits a peak for an intermediate value of the quenching 
rate. We show that though discord and also the mutual information exhibit a similar behavior with 
respect to the quenching rate to that of concurrence or negativity following an identical quenching, 
there are quantitative differences. Our studies indicate that like concurrence, discord also exhibits 
a power law scaling with the rate of quenching in the limit of slow quenching though it may not 
be expressible in a closed power law form. We also explore the behavior of discord on quenching 
linearly across a quantum multicritical point (MCP) and observe a scaling similar to that of the 
defect density. 
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I. INTRODUCTION 

Quantum Phase Transitions(QPT) [H-Q driven by 
quantum fluctuations arising due to the change of a pa- 
rameter in the Hamiltonian at absolute zero temperature 
have been studied extensively. A QPT is characterised 
by a fundamental change in the symmetry of the ground 
state of a quantum many-body system and is also associ- 
ated with diverging correlation length as well as a diverg- 
ing relaxation time at the Quantum Critical Point(QCP). 
Over last few years, numerous efforts have been directed 
to understanding the connection between Quantum In- 
formation and QPTs In recent years, QPTs have 
been observed experimentally in a large number of sys- 
tems, for example in optical lattices where a Mott insu- 
lator to superfluid transition is observed 



The entanglement between two spins is a measure of 
the correlations between them Q and is usually quanti- 
fied in terms of quantities like concurrence and negativity 
[isj . For a transverse field Ising model concurrence 
has been found to maximize close to the QCP and its 
derivatives show scaling behavior characteristics of that 
QCP 0. However, a different and significant measure 
other than the entanglement, namely the "Quantum Dis- 
cord" was introduced by Olliver and Zurek jl4| which ex- 
ploits the fact that different quantum analogs of equiva- 
lent classical expressions can be obtained because of the 
fact that a measurement perturbs a quantum system. 
This property enables us to probe the "quantumness" of 
a system. Quantum discord which ideally is a subject of 
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interest in quantum information theory [15H20l | has been 
studied for spin systems and also close to QCPs [13, [2TI - 
[25I and thereby establishes a natural connection between 
these two fields. Very recently an experimental study to 
measure quantum discord using an NMR set up has been 
reported |26t . 

In this paper, we study the quantum correlations 
present in the final state of a one dimensional transverse 
XY model after quenching the system through an Ising 
critical point between two spins separated by a lattice 
spacing n and quantify it in terms of quantum discord. 
In the process, we also investigate the classical correla- 
tions and the mutual information bewteen two spins in 
the final state and study their behavior as a function 
of the rate of quenching. We compare our observations 
with the behaviour of two spin entanglement in the final 
quantum state following a similar quench as reported in 
a recent study [13] ■ We also investigate quantum dis- 
cord following a quantum quench across a multicritical 
point (MCP) along a linear path. We note that similar 
quenching studies have been carried out to establish the 
universal scaling relation of the defect density namely the 
Kibble-Zurek scaling [2^ [2^ generated following critical 
[sol - lsH quenches. In reference [131, it was established 
that concurrence does also follow the same Kibble-Zurek 
scaling relation as the defect density; this prediction was 
found to hold good also for quenching through a MCP in 
a later study [3^ . We attempt to address the same ques- 
tion related to the scaling of discord; our studies indicate 
a similar scaling. 

The organization of the rest of the paper is as follows. 
In Sec. II, we quantify discord in terms of classical corre- 
lations and quantum mutual information. In Sec. Ill, we 
compute quantum discord of transverse XY model after 
quenching the system through critical and multicritical 
points. This is followed by a discussion on our main re- 
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suits in Sec. IV. We present some concluding remarks in 
Sec. V. 

II. QUANTUM DISCORD 

Let us consider a classical bipartite system comprising 
of two subsystems A and B. The information associated 
with the system is quantified in terms of Shannon en- 
tropy H(j>) where p is the probability distribution of the 
system. The classical mutual information is defined as 



I{p) = H{p^) + H{p^) - H{p), 



(1) 



where H{p^),i ~ A,B stand for the entropy associated 
with the subsystem i; this can alternatively be expressed 
as 



J{p) = H{p^) - Hip\p^), 



(2) 



where H{p\p^) = H{p) — H{p^) is the conditional en- 
tropy. In the quantum context, the classical Shannon 
entropy functional gets replaced by the quantum von- 
Neumann entropy expressed in terms of the density ma- 
trix 'p' acting on the composite Hilbert space. The natural 
quantum extension of Eq. ([IJ is given by 



I{p) = + s{p^) sip). 



(3) 



The conditional entropy based on local measurement 
however alters the system. The measurement are of von 
Neumann type having a set of one dimensional projec- 
tors {Bk} that sum up to identity. Following a local 
measurement only on the subsystem B, the final state pk 
of the composite system, which is the generalization of 
the classical conditional probability, is given by 



Pk = —{i(E)Bk)p(i(E)Bk), 
Pk 



(4) 



with the probability pk = tr(/ ® Bk)p{i ® Bk) where 
/ is the identity operator for the subsystem A. Quan- 
tum conditional entropy can be defined as s{p\{Bk\) = 
J2kPk^iPk)j such that the measurement based quan- 
tum mutual information takes the form J{p\{Bk}) = 
s{p^) — s(p|{i3fc}). This expression maximized based on 
the local measurement gives the classical correlation f4ll| . 
Hence we have 



C{p) 



max 



J{p\{Bk\). 



(5) 



This line of arguments provides us with two quantum 
analogs of the classical mutual information: the origi- 
nal quantum mutual information /(p) (Eq.Q) and the 
measurement induced classical correlation (Eq.®). As 
introduced by Olliver and Zurek [ij], the difference be- 
tween these two, i.e. 



Q{p) = I{p) - C{p) 



is the quantum discord which measures the amount of 
quantumncss in the state. It is noteworthy that / rep- 
resents the total information (correlation) whereas C is 



the information gained about A as a result of a measure- 
ment on i?. If Q = 0, we conclude that the measurement 
has extracted all the information about the correlation 
between A and B, on the other hand, a non-zero Q im- 
plies that the information can not be extracted by local 
measurement and the subsystem A gets disturbed in the 
process, a phenomena not usually expected in classical 
information theory . 



III. THE MODEL AND PAIRWISE 
CORRELATIONS 

We study pairwise correlations in a one-dimensional 
spin-1/2 XY model in a transverse field with near- 
est neighbor ferromagnetic interactions described by the 
Hamihonian El-lil 



H 
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[(1 + + (1 - ^)aWr + hdW. (7) 



where ct's are the Pauli spin matrices and the subscript 
stand for the spin direction and superscript the lattice 
index. The parameter h is the magnetic field applied in 
the transverse direction and 7 measures the anisotropy 
in the in-plane interactions; 7 = 1 refers to the trans- 
verse Ising model Q]. The model can be exactly solved 
by mapping the spins to spinless fermions via a Jordan- 
Wigner transformation [44|; the phase diagram for the 
model is shown in Fig. ([1]). 

We study the behavior of quantum discord in the final 
state after quenching the system across an Ising critical 
point following the quench scheme h{t) = t/T with t going 
from —00 to 00 The diverging relaxation time close 
to the QCPs at ft. = ±1 lead to defects in the final state. 
At i — > —00, the system is at the ground state |0) where 
all the spins are aligned in the -\-z direction. At t — > 00 
the system is in an excited state in which the probabilities 
of excitation for the mode k is given by 



Pk = cxp(— 7rr7^ sin^ k). 



(8) 



In the limit r — > 00, only the modes close to the critical 
modes (fc = or fc = tt) contribute and one gets pk = 
exp(— 7r7^fc^r). 

We further extend our study by quenching the sys- 
tem across a quantum multicritical point (MCP) by ap- 
proaching along a linear path [ssf 



h{^) = l + \^{t)\sgn{t)- 7(t) = --, 

T 



(9) 



and investigate the dependence of discord. The proba- 
bility of defect formation pk [36l - [38j 



(6) p/c = exp(— 7rr(l + cos fc)^ sin^ fc). 



(10) 



We shall now calculate various elements of two-spin 
density matrix in the final paramagnetic phase of the 
Hamiltonian(I71) for spins at the sites i and j = i + n 
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FIG. 1: (Color online) The phase diagram of one dimensional 
XY model in a transverse field. The vertical bold line at /i = 1 
denotes Ising transition from ferromagnetic phase to param- 
agnetic phase. The horizontal bold line stands for anisotropic 
phase transition between ferromagnetic phase ordering in x 
and y directions. A linear path to approach the MCP 'A' is 
shown. 



using the generic form of the density matrix given by 



The expressions for C\ and ci differ for different value of 
n; and it is presented belovir for n < 6: 



Cl 



C2 



f f(l-2/3o) 

(1 - 2/3o)2/?| - 4/?|+ 

f (l-2/3o)3- 2/31/34(1 -2/3o) 

i[/36{(l - 2/3o)2 - 4/3|)} + 4/32{/5|4 
/3| - /34(1 - 2/3o)}] X [16/31/34+ 
[ (1 - 2/3o){(l - 2/3o)2 - 8/3| - 4/3|}] 



,n = 2, 



, n = 6. 



The eigen values of the density matrix are obtained in 
terms of the correlators c,s [22,S31 as 



Ao 
Ai 

A2 

As 



l[(l + C3) + y'4cl + (ci-C2)2], 

i[(i + c3)-y'4cl + (ci-c2)2], 

i[(l - C3) + (ci + C2)], and 
i[(l-C3)-(ci+C2)], 



(16) 



which can be expressed entirely in terms of /3's using the 
equations dUl) and (fTS)) . 



030-3 (g) (T3 + C4/* 



where c\ = (ctiCT^), C2 



C5CT3 ' 



C3 



(11) 



C4 = C5 



((t|). For the Hamiltonian ([7]), the X and Y directions 
are equivalent and hence c\ = C2. The density matrix 
can also be expressed in the form 



P = 



fe? 
«S h'i 




a, 



fe'i'* 





a 



(12) 



where the matrix elements are given in terms of the two- 
spin correlation functions in the following manner: 



al = -((I±4)(l±a^+")) = 1 + C3±2c4, 



= -((l±4)(lTar"))=l-C3, 
^1(2) - 



(13) 



We note that the up-down symmetry of the Hamiltonian 
simplifies the density matrix and some of the elements 
vanish [4^ . Defining a quantity [23, [12] 



/3n 



/O 7^ 

one gets 



Pfc cos(nA:), 



C4 = C5 
C3 



1 - 2/3o 



4/3^ 



(14) 



(15) 



IV. RESULTS 

In this section, we present results of the pairwisc cor- 
relations in the final state of the spin chain as a function 
of the quenching rate t~^. Using Eq. we note that 
fin ~ for odd n as is invariant under k ^ ir — k. 
At the same time, {aj.a'^^^) = 6" = for all n since 
the expectation values of a pair of fermionic annihila- 
tion or creation operators do always vanish. Moreover, 
((t!|_(t^") = 62 = for odd n since t he q uantities 62 are 
odd under the Z2 transformation[23, [l^. On the other 
hand, 63 ci + C2 for even n. 

The variation of mutual Information /, the classical 
correlation C and the quantum discord Q = I ^ C with 
T are therefore studied for both critical and multicritical 
quenches ([9]) for even n. Let us rename the spin i as 
subsystem A and j as subsystem _B. The reduced density 
matrix for the subsystems A and B can be expressed as 

PA = ^{r (E}P +C4r (E)<ji), and 



with eigenvalues 

A4 = -^(^ + 04), and 



(17) 



A., 



2(1 -C4). 



(18) 



The total mutual information I{p) is expressed terms 
of von Neumann entropies, which when substituted in 
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Eq. (HI gives 

3 

I{p) = s{pa) + s{pb) - X! ^" ^°S2 Ao 



(19) 



Q=0 



where s{pa) = s{pb) = — A4 logj A4 — A5 log2 A5. To cal- 
culate the classical correlation, we introduce a set of pro- 
jector for local measurement on the subsystem B given 
by Bk = VnkVf where Uk = \k}{k\ : fc = +, - is 
the set of projectors on the computational basis |+) = 
^(|0) + |1)), I-) = ^(|0) - |1)) and V G U{2) where V 
is parametrized over a Bloch sphere given by 



cos : 



e„-i4> 



r,i4> 



(20) 



where the polar angle 9 lies between and tt and the 
azimuthal angle (j) is from to 2tt. Following a technique 
used in the reference [23|, we can obtain the classical 
correlation by maximizing 



= s{pa) - s(p+), 

where p+ is the density matrix for the outcome \k) 
Below we summarize the final results; e.g., for n = 
get 



ci = C2 = ^(l-2/3o), 



(21) 

= !+)■ 
2, we 



C3 



(l-2/3o)'-4/32, c4 = l-2/3o. 



(22) 



and 

C = -(l-^o)log2(l-/3o)-^olog2(/3o) 
+ 1 (l-(l-2/3o)yi + ^ I X 



-{1 - (1 - 2/3o 



l0g2 



+ - l + (l-2/3o)\/l + 



l0g2 




(25) 



Similarly one can obtain the expressions for n ~ A and 
n = 6 using the appropriate equations. We note that / 
and C and hence Q = I—C depends entirely on /3's which 
arc in turn dependent on the quench rate through 
the defect density p^. In deriving the above expressions, 
we have used 7=1, however, qualitatively the mathe- 
matical form of eigenvalues presented in Eq. (|16p remain 
unaltered though /?'s are modified for 7 7^ 1. 




ci = C2 = -[/?6{(l-2/3o)'-4/32)} 

+ 4/32{/32'+/34-/34(l-2/3o)}] 

X [16/32/34 + (l-2/3o){(l-2/3o)'- 8/32' 

C3 = (l-2/3o)'-4/32, c4 = l-2/3o. 



4/3|}], 



(23) 



The exact expressions for mutual information and clas- 
sical correlation for 71 = 2 is given below. 

/ = -2(l-/3o)log2(l-/3o) 

+ ((1 - /3o)' - 131) log2((l - - /?!) 
- 2/3o log2(/3o) + (/3o' - /?2) log2(/3o - Z?') 
+ i{4/3o(l " /3o) + 4/3| + /32(1 - 2/3o)} x 



l0g2 



-{4/3o(l-/3o)+4/32-f /32(l-2/3o)} 



+ -{4/3o(l - /3o) + 4/3| -Ml- 2/3o)} x 



l0g2 



-{4/3o(l-/3o)+4/3|-/32(l-2/3o)} 



(24) 



FIG. 2: Quantum discord Q as a function of r n = 2 (solid 
line), 4 (dotted line) and 6(dashed line) in the final state fol- 
lowing a linear quench across the Ising critical point with 
7 = 1. Inset shows the variation of concurrence (Cnc) for 
same parameter values as reported in [27j . 

Fig. 1^ shows the variation of quantum discord Q with 
r for n = 2, 4, 6 for a quench across the Ising critical 
point. As expected, Q vanishes in both the limits r — >■ 
and T — >■ cx); the final state is nearly a direct product state 
in either cases, discord initially increases with increas- 
ing r, and starts decreasing monotonically after reaching 
a peak at r = r"*. As n increases, r™ shifts towards 
the right. A similar behavior is observed for /. Fig. ([3]) 
shows that the classical correlations also exhibit a qual- 
itatively identical variation with r though it is smaller 
in magnitude in comparison to discord. Surprisingly, the 
classical correlations however show some fluctuating be- 
haviour for T — > followed by the monotonie increase 
(see inset Fig. ([3])). The value of C is found to be one or- 
der of magnitude less than / implying that correlations 



5 



0.0025 rr 




50 100 150 200 250 



T 



FIG. 3: Variation of classical correlation with r for n = 2 
(solid line) ,4 (dotted line) and 6(dashed line). Inset shows 
small peaks for r — >■ and it increases monotonically. 
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FIG. 4: The variation of discord with r for ?i = 2 and 7=1 
shown on a log-scale; the slope is « —0.5. Inset shows the 
variation of Concurrence with r [SJl which shows a similar 
scaling. 

present in the system are mainly quantum mechanical. 
Though C peaks at a larger t in comaprison to / for a 
given n, it does not substantially influence the behavior 
of Q as quantum correlations apparently dominate over 
classical correlations. 

The variation of concurrence in the final state for a 
quench across an Ising critical point with t has been 
studied recently [l^] and the comparison is shown in 
Fig. Although the variation of discord and concur- 
rence are qualitatively similar, we emphasize following 
differences. The magnitude of discord is less than that of 
concurrence for the same n by one order and it shows a 
peak at a value of r which is very small in comparison to 
the corresponding r for concurrence. We conclude that 
the measurement based approach employed in calculat- 
ing discord provides a quantitatively different result for 
quantum correlations. Moreover, the study on concur- 
rence [i^l indicates the existence of a threshold value of 
T above which the bipartite entanglement is generated. 
On the contrary, discord is non-zero for all r as wc ob- 
serve negligible shift close to t = for different n as 




-2 - 
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logt 

FIG. 5: The variation of discord for n — 2 following a multi- 
critical quench along a linear path with corresponding slope 
is « —0.17 which matches well with the value of the expo- 
nent — 1/6 obtained for the defect density [s^. In the inset, 
we show the similar variation for concurrence with r [s^ and 
the slope is « —0.13 which is in close agreement with the 
exponent —1/6. 

shown in Fig. ([5]). 

We would now like to investigate the scaling of Q or / 
(since magnitude of C is relatively smaller) as a function 
of T like concurrence which was found to scale as l/y/r 
in the limit of large r for quenching through the Ising 
critical point [l^l ■ To explore the scaling of / in the limit 
T — > (X), we analyse the asymptotic behavior of the terms 
present in Eq. While the first two terms taken to- 

gether show a faster decay as l/r, the total contribution 
from the remaining four terms scales numerically as 1 /y/r 
and hence dominates when r — >■ oo. Although, a closed 
power-law form is not obtained, our studies apparently 
points to the fact that discord does also satisfy a scaling 
analogous to concurrence or defect density. 

We further verify this claim by investigating the scaling 
of discord following a linear quenching across the MCP 
'A' (see Eq. (9)) for which the defect density scales as 
T^^/^ [36l - [38j . In Fig. dS]), we compare the slope of dis- 
cord with that of concurrence with respect to r on a 
logarithmic scale and find a good matching which again 
indicates that the scaling of discord is likely to be same 
as that of the defect density. 



V. CONCLUSION 

We have studied quantum discord and mutual infor- 
mation and their dependence on the quenching rate for 
two spins separated by n lattice sites in the final state 
of a transverse XY spin chain following a slow quench 
across a quantum critical and a multicritical point. Our 
studies show that the quantum discord and mutual infor- 
mation exhibit qualitatively similar behavior to that of 
concurrence as reported before [l^j; all these quantities 
are in fact determined in terms of the eigenvalues of the 
two-spin density matrix. However, we do also highlight 
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differences; discord is smaller in magnitude in comparison 
to concurrence and is non-zero even in the limit of t — > 
for any n unlike concurrence which sets in at a threshold 
value of T especially for large n. The classical correlation 
is found to be smaller than the mutual information by 
one order of magnitude suggesting a stronger quantum 
nature of correlations. Finally, our studies indicate that 
for a linear quenching through a MCP, discord shows a 
power-law scaling with the rate of quenching for large r 
in the similar fashion as the defect density. This observa- 



tion apparently suggests that the scaling of discord in the 
final state in some sense is given in terms of that of the 
defect density for both critical and multicritical quenches 
though a closed power law form is not obtained. 
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